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The Ising transition in 2D simplicial quantum gravity - can Regge 
calculus be right?* 
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We report a high statistics simulation of Ising spins coupled to 2D quantum gravity in the Regge calculus 
approach using triangulated tori with up to 512 2 vertices. For the constant area ensemble and the dl/l functional 
measure we definitively can exclude the critical exponents of the Ising phase transition as predicted for dynamically 
triangulated surfaces. We rather find clear evidence that the critical exponents agree with the Onsager values for 
static regular lattices, independent of the coupling strength of an R 2 interaction term. For exploratory simulations 
using the lattice version of the Misner measure the situation is less clear. 



1. INTRODUCTION 

Ising spins coupled to two-dimensional (2D) 
Euclidean quantum gravity in the dynamically 
triangulated surface (DTS) formulation exhibit 
a third-order phase transition with critical ex- 
ponents which are very different from the Onsa- 
ger values for static 2D lattices. This follows from 
exact results for matrix models [1] and numerical 
simulations [2], and is in perfect agreement with 
predictions from conformal field theory for mat- 
ter of central charge c = 1/2 in the continuum 
limit [3]. 

An alternative formulation is the Regge calcu- 
lus [4] where the connectivity of the discretized 
surface is fixed but the link lengths vary. Recent 
Monte Carlo (MC) simulations [5] of pure 2D gra- 
vity were interpreted in accordance with the DTS 
results. When coupled to Ising spins, however, 
evidence for Onsager critical exponents was clai- 
med [5]. Since this would be a severe failure of 
the Regge approach we decided to investigate this 
question again with much higher precision [6]. Fi- 
nally it should be remarked that already the con- 
sistency of pure gravity with the dl/l measure was 
recently questioned in Ref.[7]. 
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2. MODEL AND SIMULATION 

Employing the usual transcription [8] of conti- 
nuum quantities like the metric g and the scalar 
curvature R, the Regge partition function reads 

Z = J2 fvfi(l)exp(-I(l)-KE(l, S )), (1) 

W J 
where 

/(/) = WAA+a^) (2) 
is the discretized gravitational action and 

e(i, S )=\j: a >s(' !i t Il ) 2 (3) 

(a) v 13 j 

is the energy of Ising spins s 8 - = ±1, which are 
located at the N = L 2 vertices i of a triangulated 
torus with fixed coordination number q = 6. Here 
6{ = 2-7T — X2tDi ^' W i s the deficit angle and Ai = 
^2 tD i Af/3 is the baricentric area, where 6i(t) de- 
notes the dihedral angle at vertex i and At the 
area of a triangle t. Finally, Aij = ^2 t D{ij) At/3 
is the area associated with a link (ij). In the 
first set of simulations we used the commonly em- 
ployed simplest scale invariant functional measure 
T>ij,(1) = Yl(ij) dhj/hj. To investigate the depen- 
dence of our results on this choice we also per- 
formed exploratory simulations using the physi- 
cally more motivated lattice version of the Misner 
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Critical parameters of the Ising transition. 
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0.123(4) 
0.127(3) 
0.123(2) 


-0.06(13) 

-0.06(5) 

-0.07(9) 



measure T>[i(l) = Yl{ij) dhj / \/ Aij [5]. Since both 
measures are scale invariant the total area is ea- 
sily kept fixed at its initial value A = Ai = N 
by rescaling all link lengths when proposing a link 
update. In our simulations one MC step consisted 
of a single-hit Metropolis lattice sweep to update 
the link lengths lij followed by a single-cluster flip 
to update (a fraction of) the spins s 8 -. In each run 
we recorded the energy density e = E/A and the 
magnetization density m = AiSi/A in a time- 
series file. The statistical errors are computed by 
the jack-knife method. 

3. RESULTS 

3.1. dl/l measure 

Here we concentrated on finite-size scaling 
(FSS) analyses for the R 2 couplings a = 0, 
0.001, and 0.1, and performed long simulations 
with about 50 000 measurements near the critical 
coupling K c for various lattice sizes up to L = 512 
[6]. We did not encounter any equilibration pro- 
blems and the autocorrelation times of e and m 
turned out to be about 1 — 4 measurements for 
all lattice sizes. By applying standard reweigh- 
ting techniques we first determined the maxima 
of the susceptibility, \ = A((m 2 ) — (|m|) 2 ), the 
specific heat, C = K 2 A((e 2 ) — (e) 2 ), as well as 
of d{\m\) I dK , d\n{\m\)/dK , and d\n{m 2 )/dK . 
This defines five sequences of pseudo-transition 
points A' max (_L) for which the scaling variable 
x = (A' max (_L) — K^L 1 !" should approach a con- 
stant for large L. 

The critical exponent v can then be estima- 
ted from (linear) least square fits to the FSS An- 
satz dU L /dK = L^-^foix) or dln(\m\P)/dK ^ 
L 1 !" fp(x) to the data at the various A' max (_L), 
where Ul = 1 — (m 4 )/3(m 2 ) 2 . By averaging these 
estimates we obtain the values given in Table 1 
which are all compatible with the Onsager value 
of v = 1. Assuming thus v = 1 we next de- 



termined K c from fits to A' max (_L) = K c + c/L. 
Another, more precise method is to analyze the 
crossing points K x of the curves Ul(K) with L 
and L' = bL. Combining this information we ob- 
tained the values quoted in Table 1 which were 
used in further analyses such as the asymptotic 
limit U* of Ul(K c ). Within error bars our values 
in Table 1 agree with the very precise estimate 
for the regular square lattice, U* = 0.611(1) [9]. 

The exponent ratios j/u and f3/v follow from 
the FSS x = L~fl v f 3 {x) and (\m\) = L~^l v ' f 4 (x) 
or d{\m\) / dK = L^ 1- ^!" f^(x), respectively. The 
final averages are again collected in Table 1. Also 
here we see little influence of the R 2 term. While 
our estimates for j/u are slightly below the Onsa- 
ger value of 1.75, we obtain almost perfect agree- 
ment with the Onsager result f3/v = 0.125. Fi- 
nally we have checked that 3-parameter fits of the 
form C max = A + BL a l" yield values of ajv con- 
sistent with zero; compare Table 1. Similar fits of 
C at the other A' max (_L) sequences as well as at 
K c gave compatible results. 

Since the R 2 interaction term does not affect 
the critical exponents we give in Table 2 as final 
estimates the weighted average of the three si- 
mulations with different coupling a. For the dl/l 
measure we definitely can exclude the DTS ex- 
ponents and find strong evidence for the Onsager 
universality class. We conclude by mentioning 
that in the meantime we have also confirmed the 
failure of the Regge approach with dl/l measure 
already for pure gravity [10]. 
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Critical exponents of the Ising transition. 
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Figure 1. Distribution of local deficit angles. 

analyses of simulations with the Misner measure 
do not exclude the DTS exponents, we still have 
the hope that the Regge approach can survive as 
an alternative tool to describe quantum gravity. 
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